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Chern-Simons (CS) modified gravity is a four-dimensional, effective theory that descends both 
from string theory and loop quantum gravity, and that corrects the Einstein-Hilbert action by 
adding the product of a scalar field and the parity- violating, Pontryagin density. The Chern-Simons 
modification deforms the gravitational field of spinning black holes, which is now described by a 
modified Kerr geometry whose multipole moments deviate from the Kerr ones only at the fourth 
multipole, ^ = 4. This paper investigates possible signatures of this theory in the gravitational 
wave emission produced in the inspiral of stellar compact objects into massive black holes, both 
for intermediate- and extreme-mass ratios. We use the semi-relativistic approximation, where the 
trajectory of the small compact object is modeled via geodesies of the massive black hole geometry, 
while the gravitational waveforms are obtained from a multipolar decomposition of the radiative 
field. The main Chern-Simons corrections to the waveforms arise from modifications to the geodesic 
trajectories, which in turn are due to changes to the massive black hole geometry, and manifest 
themselves as an accumulating dephasing relative to the general relativistic case. We also explore 
the propagation and the stress-energy tensor of gravitational waves in this theory, using the short- 
wavelength approximation. We find that, although this tensor has the same form as in General 
Relativity, the energy and angular momentum balance laws are indeed modified through the stress- 
energy tensor of the Chern-Simons scalar field. These balance laws could be used to describe the 
inspiral through adiabatic changes in the orbital parameters, which in turn would enhance the 
dephasing effect. Gravitational- wave observations of intermediate- or extreme-mass ratio inspirals 
with advanced ground detectors or with the Laser Interferometer Space Antenna (LISA) could use 
such dephasing to test the dynamical theory to unprecedented levels, thus beginning the era of 
gravitational wave tests of effective quantum gravity theories. 
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I. INTRODUCTION 

The experimental verification of symmetry breaking is 
one of the most powerful tools to understand in which di- 
rection to extend the current canon toward more funda- 
mental physical theories. For example, the experimental 
confirmation of the violation of charge conjugation, par- 
ity transformation and time-reversal symmetries in ele- 
mentary particle interactions forced the improvement of 
the quantum field theory of particles into what is today 
the standard model. Similarly, violation of symmetries 
in gravitational interactions can push toward generaliza- 
tions of General Relativity (GR) by providing the first 
experimental evidence of high-energy extensions. 

Gravitational parity violation can be tuned by the in- 
clusion of a Pontryagin term in the Einstein-Hilbert ac- 
tion, which defines an effective, four-dimensional gravi- 
tational theory: Chern-Simons (CS) modified gravity 
In fact, the inclusion of such a term in the action 
is inescapable in four-dimensional compactifications of 
perturbative string theory (i.e. Type I, lib, Heterotic, 
etc.) due to the Green-Schwarz anomaly-canceling mech- 
anism 0. This fact can also be extended to the non- 
perturbative sector in the presence of Ramond-Ramond 
scalars (D-instanton charges) due to duality symme- 
tries 3 . Such a term also arises naturally in loop quan- 



tum gravity when the Barbero-Immirzi parameter is pro- 
moted to a scalar field coupled to the Nieh-Yan invari- 
ant 

The action of Dynamical CS Modified Gravity (DC- 
SMG) consists of the Einstein-Hilbert term plus the prod- 
uct of a scalar field and the Pontryagin density (the con- 
traction of the Riemann curvature tensor with its dual), 
plus the action for this scalar field and/or other mat- 
ter fields. The correction proportional to the Pontryagin 
density modifies the field equations for the metric com- 
ponents by adding two extra terms to the Einstein equa- 
tions: the so-called C-tensor and an stress-energy tensor 
for the scalar field. The C-tensor depends on derivatives 
of the CS scalar and the contraction of the Levi-Civita 
tensor with covariant derivatives of the Ricci tensor and 
the dual Riemann tensor. In addition, the variation of 
the action with respect to the CS scalar field leads to an 
equation of motion for this field, which is sourced by the 
Pontryagin density. 

The CS gravitational modification has been investi- 
gated mostly in the non-dynamical framework, in which 
the scalar field is non-dynamical (there is no kinetic 
term for it in the action), and hence it is assumed to 
be an a priori prescribed spacetime function. Such 
studies include an analysis of exact solutions [l|, , ap- 
proximate solutions P, 0, H, S Hi [m J matter interac- 
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tions Mill jcosmology [1, [H, [H, , and astrophysical 
tests Non-dynamical CS modified gravity has 

been shown to be theoretically problematic in relation to 
Schwarzschild black hole perturbation theory the 
existence of stationary and axisymmctric solutions [2]| . 
and the uniqueness of solutions of the theory [l^l • 

The detailed study of the dynamical formulation of 
CS modified gravity has only recently begun. This pa- 
per is the second in a series that deals with the details 
of the dynamics of CS modified, spinning black holes. 
In the first paper [22] . henceforth Paper I, an approx- 
imate solution was found for a spinning black hole, us- 
ing the slow-rotation approximation and a small-coupling 
approximation. The first type of approximation restricts 
attention to black holes with small angular momentum 
per unit mass, while the second one allows one to search 
for small CS deformations of known GR solutions. The 
new solution corresponds to a deformation of the Kerr 
metric whose deviations fall off with a high power of the 
distance to the black hole. 

In this paper, we concentrate on the study of inter- 
mediate and extreme-mass ratio inspirals (IMRIs and 
EMRIs respectively) in the context of DCSMG. Such 
systems consist of a small compact object (SCO) (with 
masses in the range 1 — 30 A/q) orbiting around a (spin- 
ning) massive black hole (MBH), (with masses in the 
range 10'* — lO^M©) in the case of EMRIs, and an 
intermediate-mass black hole (IMBH) (with masses in the 
range 10^ — 10'' Mq-, see [2^ for a review on the evidence 
of the existence of IMBHs) in the case of IMRIs. Another 
IMRI possibility would be that of an IMBH falling into a 
MBH, a system with obviously also an intermediate mass 
ratio (see j24y and references therein). The mass ratios 
involved are then in the range (10^^ — 10^^) for IMRIs 
and (10-4 - lO"'^) for EMRIs. 

EMRIs (and IMRIs involving a IMBH-MBH binary) 
are important sources of_gravitational waves (GWs) for 
future space detectors 25] as the Laser Interferometer 
Space Antenna (LISA) 26, 27, 28, 29, 30], whereas IMRIs 
are important sources for second generation ground based 
detectors (see [3l|,^[), such as Advanced LIGO Q and 
Advanced VIRGO ■ 34| . and for future planned third gen- 
eration detectors as the Einstein Telescope [35|- Both 
IMRIs and EMRIs are high-precision sources that can 
produce a high number of detectable GW cycles. Hence, 
GW observations of these systems will have a strong 
impact on astrophysics, cosmology, and fundamental 
physics m, [13, ' the latter of which we shall be con- 
cerned with here (see 's^ for a detailed account of IMRI 
and EMRI astrophysics and related science). 

The construction of IMRI/EMRI waveforms for data 
analysis purposes is a difficult problem since the accuracy 
needed for detection and extraction of physical parame- 
ters is quite high. For the case of EMRIs, techniques for 
constructing sufficiently accurate templates for detection 
is currently underway, mainly through the use of the adi- 
abatic approximation [H El, \M, El, S, IH, \M, 113, IH, 
[49} . Methods to build sufficiently accurate templates for 



parameter estimation has not yet been fully developed, 
the main difficulty being that one requires a more precise 
treatment of the self-gravity of the SCO and its impact 
on the gravitational waveform. For IMRIs, one might 
have to incorporate the finite-size effects of the SCO on 
the waveforms, since the mass ratio is not so extreme. 

In this article, we will use a simpler method to de- 
scribe EMRIs, the so-called semi-relativistic approxima- 
tion (s^l, in which the motion of the SCO is taken to be 
geodesic on the MBH spacetime (Kerr in GR and mod- 
ified Kerr in DCSMG) and the GWs are computed us- 
ing a multipolar decomposition (5l| . In the radiation 
zone, the different multipoles are fully determined by 
the trajectory of the SCO and its derivatives. Such a 
scheme therefore neglects radiation-reaction effects that 
scale with the square of the mass ratio of the system. One 
can add these effects by using different types of approx- 
imations that provide expressions for the change of the 
constants of motion in terms of the properties of the GWs 
emitted [H, [H, [H, [5^ . Although the semi-relativistic 
approximation is probably not sufficient to create a tem- 
plate bank for parameter estimation, it will suffice here 
to understand the type of corrections induced in CS mod- 
ified EMRIs and IMRIs. 

We begin by considering the trajectory of the SCO in 
the MBH background. We find that the motion of a test 
particle in DCSMG is exactly described by the geodesic 
equations, as in GR. Trajectory modifications are thus 
produced entirely by the modified geometry of the MBH 
background, which is here described by the modified Kerr 
solution found in Paper I. These modifications include 
corrections to the geodesic equations and to the funda- 
mental frequencies of motion. Since the CS correction to 
the MBH background is a high-order, strong- field effect, 
the multipoles of the hole are not corrected up to the 
£ — i multipole, suggesting that a test of the Kerr ge- 
ometry might be difficult with LISA. This theory is thus 
an interesting illustration of how the geometry of a MBH 
can be changed by physically well-motivated curvature 
corrections to the Einstcin-Hilbert action. 

We continue with a discussion of GW generation in 
DCSMG, which for EMRIs and IMRIs can be studied 
using black hole perturbation theory. We see that the 
perturbations, and hence the GW emission, are indeed 
corrected by the C-tensor and the stress-energy tensor of 
the CS scalar field. Expanding in the mass ratio and in 
the CS coupling constants that act as CS deformation 
parameters, we find that the second-order corrections to 
GW emission are the usual GR radiation-reaction effect 
(proportional to the square of the mass ratio) and a new 
CS correction (proportional to the product of the mass- 
ratio and the CS coupling constants). Depending on 
the strength of these constants, this modification could 
be much smaller, comparable, or bigger than radiation- 
reaction effects. In this paper, as a first exploratory study 
of these questions, we shall neglect these second-order ef- 
fects, and thus, we shall model GW emission through 
the usual GR multipolar expansion. Modifications to 
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GW emission, hence, originate from corrections to the 
geodesic trajectories of the SCO, which in turn arise due 
to the modified MBH geometry. 

With these tools at hand, we then proceed to solve the 
geodesic equations for a set of EMRI systems and com- 
pare the GWs generated in GR with those generated in 
DCSMG. Figure [T] shows the changes in the waveforms, 
essentially an accumulating dephasing, for the last eleven 
minutes of the plus-polarized waveform after 128 days of 
evolution for the system described in the caption. The 
GR waveforms are here denoted with a dotted blue line, 
while the DCSMG waveforms are denoted with a black 
solid line. The amount of GS dephasing depends on the 
specific type of orbit and the magnitude of the GS cor- 
rection, where we generically find a strong effect on the 
waveform for orbits that spend the longest close to the 
MBH. 
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FIG. 1; Plus-polarized waveform as a function of time using 
a Kerr background (dotted blue) and the CS deformed metric 
(solid black). The pericenter is at r = 6M (with M = 4.5 x 
10''A^q) the inclination angle is 0.3 radians, the spin angular 
momentum is J/M^ = 0.4, while the CS coupling strength 
is set at ^/M* = 0.4. For more details on this system see 
Sec. IV Al and for further details on the precise definition of 
the CS coupling strength see Sec. [Ill 

Our study suggests that a GW observation of highly 
relativistic IMRIs and EMRIs could place constraints on 
the dynamical theory that are orders of magnitude larger 
than current binary pulsar ones (2^ . This constraint de- 
pends strongly on how relativistic the system is, as well as 
on the signal-to- noise ratio of the event, and on the total 
mass of the system. Including radiation-reaction effects 
should allow us to search, for the first time, for radiative 
effects associated with the CS scalar field in the neigh- 
borhood of compact objects, which is simply not possible 
with binary pulsars (see [56l . |57[ for detailed information 
on tests of gravitational theories, in particular for tests 
with pulsars). To confirm these expectations a detailed 
analysis using appropriate data analysis tools is required, 
but shall not be performed here. 

Remarkably, a physically well-motivated curvature- 
squared correction to GR, like the one studied here. 



seems to lead to spinning black hole solutions and GWs 
that resemble the GR prediction quite closely. This 
fact suggests that GW detectors should be able to de- 
tect GWs, irrespective of whether GR or GS waveform 
templates are employed. As long as the precise nature 
of the massive compact object or the structure of the 
true gravitational theory does not modify the GR pre- 
dictions enough, as we show is the case in DCSMG, the 
current data analysis algorithms should be able to ex- 
tract the signals by using purely GR templates and per- 
forming incoherent searches of the data in short segments 
(e.g. roughly three weeks). However, in our case, the GS 
modification may prevent us from connecting these seg- 
ments together and associate them to a specific EMRI or 
IMRI system. 

Although the study of GWs emitted by EMRIs asso- 
ciated to non-Kerr backgrounds is not new [sl, [s^ , the 
analysis presented here is the first one to consider an al- 
ternative to the Kerr metric that derives from a concrete 
and physically well-motivated alternative theory of grav- 
ity (based on high-order curvature corrections to GR) . In 
contrast to arbitrary quasi-Kerr deformations, the black 
hole solution considered here deviates from Kerr signif- 
icantly only in the strong field region and only through 
gravitomagnetic components. Thus, mismatches between 
GW signals of EMRIs in this background (and in DC- 
SMG) and GWs in Kerr (and in GR) are dominated ex- 
clusively by truly strong-field effects and not by sublead- 
ing, weak-field corrections. Tests of such modified theo- 
ries seem to require not only the detection of a sufficiently 
large number of segments, but also the reconstruction of 
a unique GW signal. 

We conclude this paper with a discussion of the role of 
radiation-reaction effects in DCSMG through the short- 
wavelength approximation. In particular, we investi- 
gate how one could introduce these effects in the adia- 
batic approximation (that is, when the rate of change 
of the orbital parameters due to GW emission is much 
smaller that the typical orbital periods) in order to con- 
struct truly inspiraling IMRI/EMRI waveforms in DC- 
SMG. We find that the effective GW stress-energy ten- 
sor, the Isaacson tensor, and thus the GW stress-energy 
loss, is the same in DCSMG as in GR. The backreac- 
tion of this loss on the background, however, is differ- 
ent in DCSMG because GW losses must here be supple- 
mented by the stress-energy distribution of the CS scalar 
field. Truly inspiraling IMRI/EMRI waveforms can then 
be constructed in DCSMG by supplementing the semi- 
relativistic approximation with adiabatic changes of the 
orbital parameters due to both GW and CS scalar field 
stress-energy losses. 

The details of this study are organized as follows: 
Sec. [H] discusses the basics of DCSMG, its linearized ver- 
sion, the study of GW polarization in this theory, and 
reviews the small-coupling approximation; Sec. IIIII de- 
scribes test-particle motion in DCSMG, including the ex- 
pression of the modified Kerr solution, and the derivation 
of the equations for timelike geodesies. We also discuss 
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how the fundamental frequencies of geodesic motion are 
modified due to the changes in the MBH metric; Sec. IIVI 
describes GW generation in the modified theory and the 
type of semi-relativistic approximation that we use; in 
Sec.|V]we apply the previous results to construct numer- 
ically the trajectories and waveforms; Sec. IVII discusses 
GW propagation in the short-wavelength approximation 
and shows that the GW energy momentum tensor and 
emission in DCSMG is the same as in GR. We also cal- 
culate the stress-energy emission due to the CS scalar 
field and the modification to the background (averaged) 
scalar field due to terms quadratic in the perturbation 
of the Riemann tensor; Sec. IVIII concludes and discusses 
future research directions. 

We use the following conventions throughout this work. 
Greek letters and a semicolon are used to denote in- 
dices and covariant differentiation respectively on the 4- 
dimensional spacetime. We denote covariant differentia- 
tion with respect to the background metric by V^i?^ or 
by B^^^. Partial differentiation with respect to the coor- 
dinate is denoted as d^B^ or i?^^. Symmetrization 
and antisymmetrization is denoted with parenthesis and 
square brackets around the indices respectively, such as 
\^.) ■■= + ^.m]/2 and ^[^^] := [A^, - A,^]/2. We 
use the metric signature (— , and geometric units 

in which G = c = 1. 



II. CHERN-SIMONS MODIFIED GRAVITY 

In this section we describe the formulation of CS modi- 
fied gravity that we shall employ and establish some basic 
notation (see also Paper I). We begin with a discussion of 
the basics of the modified theory, but refer the reader to 
Paper I or the upcoming review [iOl for further details. 
We continue with a review of the linearized theory in 
DCSMG. Then we study basic properties of GW polar- 
ization in DCSMG, illustrated with plane waves. Finally 
we introduce the small-coupling approximation. 



A. General Formulation 

The starting point is the action of DCSMG: 



(1) 



where 



(fx /=g § *RR, 




(2) 



where k = (IGttG)"^ is the gravitational constant; a is 
the coupling constant of the CS scalar field i9 with the 
parity- violating Pontryagin density, *R i?, given by 

RR — Rap-iS ^ ~ 2 ^ap-yS^ ^iy ^ I'^j 

where the asterisk denotes the dual tensor, constructed 
using the antisymmetric Levi-Civita tensor g"/''^''; j3 is 
a constant that determines the gravitational strength of 
the CS scalar field stress-energy distribution. The dif- 
ferent terms in Eq. ^ correspond to the following: the 
first one is the Einstein-Hilbert action; the second one is 
the CS gravitational correction; the third one is the CS 
scalar field action, which contains a kinetic and a poten- 
tial term V{d), both of which distinguishes the dynamical 
formulation from previous ones; and the fourth one is the 
action corresponding to matter degrees of freedom. 

Upon variation of the action with respect to the met- 
ric and the CS scalar, we obtain the field equations of 
DCSMG: 



-C,.„ - 



1 

2k 



T, 



l_Ll' 



f3ad = (3 



dV 

M 



a 



*RR, 



(4) 
(5) 



where T^'J"' is the matter stress-energy tensor and T^? is 
the stress-energy of the CS scalar field, given by 



(V^i9)(V,i?) 



(6) 

In this work we will study only the case in which the 
potential term associated with the CS scalar field is as- 
sumed to vanish^ {V = 0). The tensor (the so- 
called C-tensor) in Eq. ([4]) can be split into two parts, 
Ct"" = Cf" + C^", where 



5 > 



c: 



a(3 



(7) 



As we have mentioned, there are two conceptually dis- 
tinct CS modified gravity theories: a dynamical version 
and a non-dynamical version. In the former, the quan- 
tities a and are arbitrary, and the field equations are 
the ones we have given in Eqs. Q and (jS]). The non- 
dynamical version is characterized by the choice (3 = 
at the level of the action, and thus the evolution equa- 
tion for the CS scalar becomes a differential constraint on 
the space of allowed solutions, the so-called Pontryagin 
constraint *RR = 0. 



^ In string theory, the i9 field corresponds to the axion, which pos- 
sess a shift symmetry and a vanishing potential along certain flat 
directions. The axion is a moduli field, which, when stabilized 
by supersymmetry breaking |61| . acquires a non-flat potential. 
Such potential would arise at the scale of supersymmetry break- 
ing, which is much larger than the gravitational scale. 
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The non-dynamical theory presents a certain number 
of difficulties that make it less physically interesting than 
the dynamical formulation (20l.l2ll.[2^. In the former, one 
can show that single-parity perturbation modes lead to 
an overconstrained system of perturbed equations, disal- 
lowing generic quasinormal ringing [20l |. Moreover, spin- 
ning black hole solutions were found to be strongly re- 
stricted by the Pontryagin constraint, in some cases also 
leading to an overconstrained system [l^l- Finally, the 
freedom in the choice of the CS scalar field was shown to 
lead to non-unique solutions with infinite energy scalars 
in Paper I 

Most of the difficulties and arbitrariness of the non- 
dynamical theory can be avoided by adopting the dy- 
namical framework, i.e. DCSMG theory, where the free- 
dom regarding the CS scalar field reduces to the choice of 
initial conditions. The stationary state of the CS scalar 
has been seen to be independent of the initial conditions, 
being determined only by the metric through the source 
term in the evolution equation for i!} 22]. For these rea- 
sons, we choose here to study the dynamical formulation. 



B. The Linearized Theory 

In perturbation theory we can split the spacetime 
metric as the sum of a background metric g^^ (we use 
the overbar to denote objects associated with the back- 
ground) and a metric perturbation /i^^: 



and Ricci tensors are j62l |: 

*^5a/37 ^ ^S^,^af^Jf3M^^ (12) 

respectively, where — rj^'^df^di, is the flat-space 
D'Alembertian and e^i^a/s jg flat -space Levi-Civita 
symbol. The linearized C-tensor in this gauge is then 

- i(9.^^?)e-*«(^| [/i5^,|,)5-/i.)c,7A-] . (14) 

Combining these linearized expressions, we find that the 
linearized modified field equations in trace-reversed form 
become 

--t;, = □^/^^. + -(a,i9)e-'^^«(^n^/i,)5,« 



+ -(a,^79)e-^«(^ [/i^^^|,)5~/i,)^,^,J ,(15) 



(16) 



(8) 



where T^^ = T^^ - (1/2)77^^T'' is the trace-reversed 
stress-energy tensor of the CS scalar-field. We see clearly 
that the metric perturbation is sourced by T^^^ which im- 
plies that the contribution of this source must also be of 
order ©(e), while the scalar field is sourced by the Pon- 
tryagin density that is of 0{t^) in the radiation zone. 
Note in this regard, that we have not yet introduced any 
perturbative expansion for the CS scalar field. 



where e is a book-keeping perturbative parameter, associ- 
ated to the smallness of the metric perturbation relative 
to the background, which we use to label the order of the 
approximation. In this section and in other parts of this 
paper, we simplify equations by using the Lorenz gauge 

l^h'^- -^^hr""^ ^ ^0, (9) 

where h — g^'^h^^ is the trace of the metric perturba- 
tion with respect to the background. Since we shall be 
at first interested in computing the DCSMG corrections 
to gravitational radiation measured by observers in the 
radiation zone, we shall restrict ourselves to a fiat back- 
ground, g^jj = T]^^, which is a good approximation suffi- 
ciently far away from the source. In the radiation zone, 
we can further restrict the gauge by supplementing it 
with the trace-free condition 

/i = 0, (10) 

leading to a transverse-traceless gauge far away from the 
source. 

In this gauge, the linearized Riemann, dual Riemann, 



C. Gravitational- Wave Polarizations 

Since we are dealing with a different theory of grav- 
ity, one must be specially careful of not assuming GW 
properties that hold only in GR. In this sense, a ques- 
tion of particular relevance that arises with alternative 
theories of gravity refers to the number of independent 
GW polarizations. DCSMG possess an extra degree of 
freedom described by the CS scalar field, and thus, ad- 
ditional GW polarizations could in principle be present. 
This issue can be investigated by following the pioneer- 
ing work of [H, [g^l, where a formalism is presented to 
study far-field gravitational radiation in any metric the- 
ories of gravity. One then finds that only six possible 
independent GW polarizations are possible, leading to a 
suitable classification (the so-called E{2) classification) 
of alternative theories. 

In this framework one focuses on the propagation of 
monochromatic plane waves in the radiation zone (i.e. 
assuming weak-fields) and their effect on test masses. 
This can be done by considering the geodesic deviation 
equation (see e.g. [S^Hl]), which describes the relative 
acceleration between nearby geodesies: 

piVV^X^") A^" = Rf^p^WUPX" , (17) 
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where is the test mass four-velocity and X^^ is the 
so-called deviation vector, which describes the displace- 
ment between test masses. Different GW polarizations 
will induce a different deviation effect on test masses, 
that can be classified via the structure of the Riemann 
tensor. Thus, the E{2) classification sorts metric theories 
based on the vanishing or non- vanishing of the Newman- 
Penrose scalars (in an appropriate Newman-Penrose null 
basis), associated with the Riemann curvature tensor, 
which can be shown to describe GW polarizations. 

Applying this procedure to DCSMG (one can use the 
formulae of the linear theory presented above), we find 
that only the two GW polarizations also present in GR 
are observable far away from sources. Such a study had 
not yet been performed in DCSMG, although several 
investigations exist in the non-dynamical theory: [I| 
proved the above result in the non-dynamical theory for 
a linearly time-dependent scalar field, while [l3| ex- 
tended this result to generic time-dependent scalars. 

To illustrate this fact in the non-linear regime, let us 
consider the propagation of plane-fronted gravitational 
waves (pp-waves) in DCSMG (see plj for pp-waves in 
the non-dynamical formalism). The line element for pp- 
waves can be written as (see e.g. [gtI): 

ds^ ^ -2dudv - H{u, x, y)du^ + dx^ -|- dy^ , (18) 

where u — t — z and v — t + z are retarded and advance 
null coordinates respectively (the waves propagate along 
the z axis), the wave fronts are given by the null planes 
u = const., and all physical information is encoded in the 
scalar H . Indeed, using a Newman-Penrose null basis [67, 
[68| adapted to pp-waves {k = , £ = d.^^ — {H/2)d^, m = 
{d^ + i dy)/ V2) one finds that the only two non- vanishing 
Newman-Penrose complex curvature scalars are: 

^,^li^l^^l + 2^^ly)H, <f,, = \idl + dl)H. (19) 

Considering an arbitrary CS scalar field, 'd = i){u, v, x, y) 
and the metric in Eq. ((T8)) . the modified field equations 
[Eqs. (HI) and ^] force us to 

d = d{u) , {dl + dl)H - ^ (20) 

where 'd{u) is an arbitrary function and = dd / du. From 
Eqs. (HH) and we see that in GR (/3 = 0) the only 
non- vanishing Weyl complex scalar is \E'4, which encodes 
the two GR GW polarizations. In DCSMG, there is an 
extra degree of freedom due to the real Ricci scalar $22 , 
which describes a transverse GW whose effect on a ring 
of test particles is to either expand or contract it main- 
taining a circular shape (a so-called breathing mode) [63j . 
The non-vanishing of ^'4 and $22 suggests naively that 
DCSMG is of type (like the Dicke-Brans- Jordan the- 
ory [69|]), but on closer inspection <I>22 is of 0{6'^), and 
thus the theory reduces to type N2 (the same as GR) 
in the weak-amplitude approximation that is required in 
the E{2) classification. 



This new polarization state can be studied more 
cleanly by considering plane waves [T^I , a special class of 
pp-waves characterized by a scalar H that is quadratic in 
X and y (^'4 and $22 independent of x and y and the 
metric has a five-dimensional Lie group of Killing sym- 
metries (67|). Terms linear in x and y can be removed by 
coordinate transformations that leave invariant the line 
element in Eq. (fT5)) . The solution to Eq. can then 
be written as 

H{u, X, y) = A[u)[x'' ~ y^) + 2B{u)xy + C[u)[x'' + y^) , 

(21) 

where A{u) and B{u) are arbitrary functions and C{u) = 
(/3/(4K))zj2. From ^ one finds that ^4 = A + iB 
and (f>22 — C*, clearly showing that the extra non-GR 
polarization corresponds to a breathing mode. In the 
weak-amplitude regime, however, iS'^ ~ 0, and the effect 
of this extra polarization is negligible, thus establishing 
the result that DCSMG only dominantly excites the same 
two GW polarization states excited in GR. 

D. The Small-Coupling Approximation 

This approximation consists of expanding the modi- 
fied field equations in the dimensionless parameter ^ = 
0(^/M'*), where f := /{/3k) and M is a characteristic 
mass associated with the particular system under con- 
sideration. To that end, we assume that C is a small 
perturbation parameter associated with the CS gravita- 
tional modifications and the physical system under study. 
As we did in the linearized theory, ^ will be also used as 
a book-keeping parameter for labeling the different per- 
turbative orders, but we will set it to unity at the end 
of our calculations. Combining this approximation with 
the linear one we can set up a two-parameter perturbative 
scheme, in ( and e, that in the context of the DCSMG 
was introduced in Paper I and we refer to this work for 
details (for a general discussion about multi-parameter 
perturbation theory see fn\ . and in the context of GR 
see [7l[7l,[7l[7l[7i[). As in Paper I, the order counting 
shall assume that /3 = 0{a) such that a/(3 — 0{1) and 
^ = 0{a/K) without loss of generality. 

The metric perturbations /i^^ and the CS scalar field 
9 can then be expanded as 

- Y.e'^C' h<^;/^ , = ^ 6"^ C' ^'•""■'^ , (22) 

a,b a,h 

where ft.j^t''''' and i?'"-'') stand for the perturbations of 
0(^°, C}') and a + 6 > 1. Introducing these expansions in 
the field equations one can then set up a boot-strapping 
method that consists of first solving the evolution equa- 
tion for the CS scalar field and then using this solution in 
the modified Einstein field equations to solve for the CS 
correction to the metric perturbations. In order to pre- 
vent the existence of metric perturbations whose domi- 
nant behavior is dramatically different from GR, we shall 
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immediately set h 



(0,a) 



for all a, as we are here search- 



ing for CS- deformations of GR solutions. 



III. TEST PARTICLE MOTION 

The first approximation to the dynamics of 
IMRI/EMRIs is to treat the SCO as a point parti- 
cle which, according to GR follows geodesic motion. In 
this section we study how this picture is modified in 
DCSMG, from the change in the geometry of the MBH 
to the change in the structure of the geodesic equations. 



A. MBH Geometry 

In GR, the geometry of a spinning MBH is described 
by the Kerr metric, but in the study of CS gravitational 
modifications it is well-known that this is no longer the 
case [2l|. Recently, the corrections to the gravitational 
field of a spinning black hole have been found [23| in DC- 
SMG using the slow- rotation approximation^, a/M ^ 1 
(where M refers to the MBH mass and a is the MBH 
spin parameter), and the small-coupling approximation 
up to second order in a/M and ^. The form of the non- 
vanishing metric components using coordinates in which 
the GR part of the metric is Kerr in Boyer-Lindquist 
coordinates {t, r, 9, <j)) is: 



1 



2Mr\ 



2Mar 
P 

5 ^ a 



■ sm 



M4 M 



12M 27AP 



7r 



El 
A 



10r2 



sin^ 6 , 



P 



(23) 



+ cos^ 0, A = r^/ 



/ = 1 - 2M/r 



where p 

and — {r^ + a^Y ~ sin^ 9. The polar angle 9 is not 
to be confused with the CS scalar field d. Technically, the 
expressions in Eq. (I23p are valid only up to second order 
m. a/M and but we have here chosen to resum the GR 
sector, by adding the appropriate high-order uncontrolled 
remainders. In addition, the expression for the CS scalar 
field (using the same approximations) is: 



- 5 a a cos(6') 



2M 
r 



,t2 



IM'I 



(24) 



that this CS scalar decays as r ^ in the far-field and thus 
it possess a finite energy. 

One could in principle attempt to calculate the 
0[{a/M)^X] corrections to the metric. These correc- 
tions, however, are difficult to find because they involve 
modifications to all metric components, since the modi- 
fied field equations must include the scalar stress-energy 
tensor to this order. We shall not consider these correc- 
tions here and work only to leading order with the metric 
of Eq. . Needless to say, if the equivalent metric for 
rapidly rotating MBHs were found, one could repeat the 
analysis of this paper for that background spacetime. 

The non-Kerrness of Eq. (|23p allows us to study DC- 
SMG modifications to the multipolar structure of a spin- 
ning black hole, as this may illustrate generic corrections 
of higher-order curvature extensions of GR. The multi- 
polar structure of Kerr is fully determined by only two 
multipole moments: the mass monopole and the current 
dipole. All others are related to these via the simple re- 
lation [78]: Mi + iSi — M {iaY , where {M^j^^Q ^ and 
{^i}e=o CO mass and current multipole moments 

respectively. In DCSMG, the leading-order modification 
to this relation occurs for the multipole, as one can 
see by employing the multipolar formalism of 51] (see 
also [79|), and noting that the only metric component 
that is CS modified scales as r^'* for M/r <C 1. 

DCSMG, therefore, preserves the idea of the no-hair 
(or two-hair) conjecture that spinning black holes are 
determined by the mass and spin parameters, but it in- 
troduces a 4-pole (or hexadecapole) correction. The rea- 
son why DCSMG does not violate the conjecture is be- 
cause, apart from mass and spin, the other parameters 
appearing in the MBH metric are introduced through 
the parameter ^, which consists of fundamental coupling 
constants of the theory that are fixed and non-tunable 
(i.e. these constants are analogous to the Newtonian grav- 
itational constant G). Instead, the modifications intro- 
duced by DCSMG are such that the standard GR formu- 
lae to relate mass and angular momentum to all multipole 
moments of the solution does not hold. 

The high ^-number of the modification suggests that a 
GW test of the MBH spacetime geometry (or test of the 
GR Kerr solution) would require a high accuracy, chal- 
lenging the abilities of present and future planned detec- 
tors. For example, LISA observations are expected to be 
able to produce accurate measurements of the first 3-5 
multipole moments of the MBH (see [13, HI Isl] for dis- 
cussions regarding this problem), which is at the bound- 
ary of DCSMG detectability. 



The metric in Eq. (|23p is stationary and axisymmetric, 
and the CS scalar field has the same symmetries. Notice 



Motion of Massive Particles 



^ After the publication of this work, other researchers [TTI have 
arrived at the same solution as that presented first in Paper I, 
suggesting a certain robustness and uniqueness of Eq. I I23I I . 



Let us now consider the equations of motion for mas- 
sive test particles in DCSMG. The starting point is 
the action of such a particle moving along a worldline 
— 2^ (A), where A parameterizes the trajectory. This 
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action is given by (see e.g. [83 



(25) 



where m is the particle mass and i'^ — dz^^/dX is the tan- 
gent to the worldhne 7. The contribution to the matter 
stress-energy tensor can be obtained by variation of S^_^.^^ 
with respect to the metric and yields 

T^£(a;'^) = m / -^u"M'^5('')[a; - z{t)] , (26) 



where t denotes proper time (which is related to A via 
dr — d\ ^ ~g^p{z)z°'zl^), = dz^ / dr is the particle 
four- velocity {g^^u^u'^ — —1), g denotes the metric de- 
terminant, and 6^'*' is the four-dimensional Dirac density 
(/dSv^6(4)(a;) = 1). 

The divergence of is 



'2 ,Q 



dr d^z 
\f^-g dT^ 



■b^^\x - z{t)] . 



(27) 



and, substituting this into the divergence of the field 
equations [Eq. (|4])], we obtain 

Since the left-hand side of this equation vanishes by 
virtue of the evolution equation [Eq. (O], all point-like 
particles must follow geodesies: 



dr^ 



0. 



(29) 



Notice that this result is completely generic and does not 
rely on any approximation scheme. 

Let us now consider the structure of the timelike 
geodesies of the spacetime described by the metric in 
Eq. (|23p . Since this metric is a small deformation of the 
Kerr solution, we can follow the same steps that lead 
to the derivation of the Kerr geodesic equations in GR 
(see e.g. |8^]) in a form suitable for numerical imple- 
mentation. The new metric is still stationary and ax- 
isymmetric, which means it possesses a timelike Killing 
vector with components = [1,0,0,0] and a spacelike 
Killing vector with components — [0,0,0,1] respec- 
tively. These two Killing vectors commute as in the Kerr 
case, and we can use them to define two conserved quan- 
tities: the energy, E = —t^Ua, and the angular momen- 
tum, L = (jf^Ua- In Boyer-Lindquist coordinates, the par- 
ticle four- velocity has components u°' = [i, r, 0, 0], where 
the overhead dots now stand for differentiation with re- 
spect to proper time r. We shall work here in reduced 
variables, where E and L are the energy and angular 
momentum per unit mass m. 

We can now use the energy and angular momentum 
definitions to solve for t and 0. To second order in the 
slow-rotation approximation, we have 



(30) 



where and (pj^ denote the corresponding expressions 
for Kerr: 



p'tK = 



—a (yttEsa 
- iaE- 



-L) 



sm 



aP 



(31) 
(32) 



and where P — E{r'^ + a^) — aL, while 6g'^^ denotes the 
CS correction: 



(3k 112r8/ 



(TOr^ + 120rM-f 189M2) 



(33) 



As before, we have here resumed the Kerr part of the 
equations, so that when we take a ^ we recover the 
exact equations for the Kerr spacetime for all a, although 
the expressions presented in this paper are only formally 
valid up to second order in a/M and C- 

Following closely the Kerr case, we look for a Killing 
tensor to construct an additional constant of motion: the 
Carter constant. We find that such a tensor can be writ- 
ten in the same form as in the Kerr case, namely 



iaf3 — ^ ^{J, 



ial3 ■ 



(34) 



where /c" and P are two null vectors given by 



A 

J2. 



.-1,0, 



1,0, 



-SgT 



(35) 



that is, they are modifications of the Kerr principal null 
directions . We have checked that they are also principal 
directions of the new metric up to the order considered 
here. One can show by direct evaluation that this tensor 
satisfies the tensor Killing equations, V^^^^,^^-, = 0, on 
the CS modified Kerr background if and only if Sg'^'^ is 
given by Eq. ([55]) . In fact, the tensor Killing equations 
allow us to add a term of the form C/r'^ to Eq. p3|) . but 
the orthogonality relation k"la = forces C — 0. 

We can now define the Carter constant in the same 
way as in the Kerr case 



(36) 



Q = i^au'^u^ - (£ - aEf 



Using this relation, together with Eqs. (|3ip and 
and the normalization condition for timelike geodesies 



-1, we can solve for f and 6. We find that 



6q/3 

the result can be written as 



rl + 2ELf5gf , 6^ = 



52 _ q2 

K , 



(37) 



where the part that corresponds to Kerr (without ex- 
panding in a/M) is 



P'rl = [{r' + a^)E-aL] 
- A 



Q+{aE- Lf + 



P 



Q - cot^ OL^ - cos^ ( 



l^E' 



(38) 
(39) 
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To summarize, Eqs. ((30)) and l|37|) provide a set of four 
first-order (in derivatives of r) decoupled geodesic equa- 
tions for the CS modified Kerr background. We note 
that these geodesic equations contain CS corrections ex- 
cept for the 9 equation, and that these corrections have 
essentially the same structure, in the sense that all of 



where the periods, denoted by A, are given by 



them are proportional to Sg'r^ [Eq. 



We remind the 



reader once more that the polar angle 9 and its evolution 
equation 9 is not to be confused with the CS scalar field 
■d and its time derivative -d. 

One can easily show that the location of the innermost 
stable circular orbit (ISCO) for equatorial orbits in the 
CS modified metric is given by [22| 



4^/6a 



^ 77V6a 



18M 5184 (3kM^ 



(40) 



where the upper (lower) signs correspond to co-rotating 
(counter-rotating) geodesies. Then, it appears that the 
CS correction works against the spin effects, which sug- 
gests that similar behavior will be present in the solution 
of the CS modified geodesic equations. 



C. Frequencies of the Orbital Motion 

The Kerr metric allows bound and stable geodesic tra- 
jectories (orbits) to be associated or decomposed in terms 
of three fundamental frequencies (see e.g. [H, [s^): 0.^ 
characterizes the radial motion (from periapsis to apoap- 
sis and back); fig characterizes the motion in the polar 
direction; and characterizes azimuthal motion. These 
frequencies are important because precessional orbital ef- 
fects are due to mismatches between them and because 
they can be used to decompose, among other things, the 
gravitational waveform in a Fourier expansion. 

Expressions for these frequencies in terms of quadra- 
tures have been obtained for Kerr in [85| . and recently 
also in (86| . where a new time coordinate is employed 
(see [13]): p^d/dr = d/d\. In terms of the time A, the 
modified geodesic equations become 

^ = r^(r,0)+ros(r), ('^y=i?,,(r)+i?cs(r). 



d(t) , ^ , , fd9 



©kW, (41) 



where (Tj^, Rj^,Qj^,^j^) are given by the right-hand 
sides of Eqs. and ^ respectively, while 

(Tj^s, i?(,3, $pg) are quantities proportional to Sg'^^. Note 
that the CS corrections only depend on r, which is a 
consequence of the linearization in a/M . 

The fundamental frequencies of the radial and polar 
motions associated with the A time are 



27r 



27r 

aF 



1 



A^ 
A5 



2lT 

a: 



27r 

aF 



(42) 




where r^^^ and r^^^j are the apocenter and pericenter val- 
ues of r respectively, and 9_^^_^ determines the interval in 
which 9 oscillates, i.e. (0„i„,7r — 6'„,i„). The frequency 
associated with the azimuthal motion is given by [86[ 



(27r)2 



2tt n2TT 

dw'' / dw'^ $[r(w'^),6'(«;^)], (45) 
Jo 



and where w^'^ = y^gX are the associated angle vari- 
ables. Therefore, only and change with respect to 
GR. 

The expressions for the fundamental frequencies with 
respect to the coordinate time t, g ^ are more involved 
since they are given by fl^gj, ^yro <h/^^ and T is [sil 



1 



dw'' / dw'^ T[r{w''),9{w^)]. (46) 



(27r)2 J, 

Inverting this relation cannot be easily done in close form, 
without previously solving the integrals numerically (see 
eg. [s^ for an analytical inversion in terms of incomplete 
elliptical integrals). We immediately see, however, that 
since T is CS corrected, so will be F, leading to CS correc- 
tions in all the frequencies associated with the coordinate 
time t. 

Since the three physical fundamental frequencies of the 
Kerr spacetime are modified in CS modified gravity, given 
a detection of quasinormal modes from a Kerr-like object, 
the inference of system parameter could lead to a system- 
atic error due to the bias of assuming GR is the correct 
underlying theory [s^l . Conversely, the detection of such 
modes with LISA may yield an interesting constraint on 
dynamical CS modified gravity. 



IV. GRAVITATIONAL WAVE GENERATION 
AND THE SEMI-RELATIVISTIC 
APPROXIMATION 

In this section we look at the structure of the modified 
field equations of DCSMG to study how gravitational- 
wave generation formulae change with respect to the GR 
ones. We also describe the semi-relativistic approxima- 
tion for the description of IMRI/EMRI systems and the 
details of the implementation that we use in this paper. 



A. GW generation in DCSMG 

Taking into account the extreme mass ratios involved, 
the systems under considerations can be well described 
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using perturbation theory around the background quan- 
tities given in Sec. IIII Al In that section, we showed that 
these background quantities contain a GR contribution 
plus perturbations in the CS couphng parameter C and 
in the black hole rotation parameter a/M. Similarly, 
here we shall expand the perturbations generated by the 
SCO (treated here as a particle) orbiting the MBH also 
in powers of a third parameter: the mass ratio of the 
system /i = m/M . 

Neglecting the slow-rotation approximation for the 
time being, the metric and the CS scalar field can be 
expanded as follows: 

n—l.rn—0 

= g^. + f^hlY^+fiChlY^ + ... , (47) 



1) = d+ E ^"C™'?^"'"^ 

n— l,m— 

= i9 + + K'?^^'^^ + ■ • ■ . (48) 

where (g^i,, ^) are given in Eqs. (|23|l and ([24|) and /i^"'™^ 
and now stand for the perturbations of C(/i", C™) 

(alternatively, one can think of /i as playing the role of e 
in Sec. Ed]). 

The equation for the metric perturbation h'j}^^'^ , in the 
Lorenz gauge [Eq. is: 

^l{ahlY^ + 2Rp/X'f^} = - It-- , (49) 

where there are no contributions from the CS scalar field 
at this order of approximation, and thus, this is a purely 
GR GW generation equation. Similarly, the equation for 

^(1,0) ig. 

(50) 

where '^TJ^^ and '^'(*i?i?) denote the first-order pertur- 
bations of the Christoffel symbols and the Pontryagin 
density respectively, about g^^ and evaluated at h'j^J^K 
That is, 1^(1'°) is generated by the metric perturbations 
/ili"' which, in turn, are generated by the motion of the 
SCO around the MBH. The first two terms on the right- 
hand side of Eq. (|50p are proportional to the spin pa- 
rameter a/M. This is not the case for the third term, 
since this is non-zero for a Schwarzschild BH [20] , which 
implies that ■d'^^'^'^ will contain terms proportional to a 
but also terms that are a- independent. 

The next correction to the metric perturbations in the 
small coupling approximation is h']^y^\ which satisfies an 



equation with the following structure 

-\h'^^f^v'^^vjY (51) 

where '"^'C^^ is the first-order perturbation of the C- 
tensor when only the metric is perturbed [see Ap- 
pendix IXTj for explicit formulae of these quantities] . Note 
that the terms on the right-hand side of Eq. ([5T|) de- 
pend only on the background (over-head barred) quan- 
tities and the (1,0) perturbations, which are determined 
by Eqs. ([iQ ]) and ([50 1 . 

The leading-order CS correction to the GW generation 
formalism is then determined by Eq. (j5ip . so let us now 
analyze its structure. The first term in this equation is 
(2a/K)C^j,[i9(^'°\ g^g.], and it contains two pieces: 

(2a/K)(V,79(i'°))e'^*"(^V„i?'')5 , (52) 
(2a/K)(V^V5Z?(i'°))*i?^(^'')". (53) 

Since the background is Kerr plus terms of order C, 
Eq. ([5^ is actually at least of 0{(^ iia/M), and we 

can then ignore it to compute /i^l^^^. The second piece 
[Eq. (|55|) ] is actually the only one in Eq. ([?T|) that is not 
proportional to the spin parameter, making it of 0{C /i), 
since the dual to the Riemann of the Schwarzschild met- 
ric is non-vanishing. The last term in Eq. ([5T|) is at least 
quadratic in a/M, making it of 0{( fi a^/M^), since 
is linear in the spin parameter. All remaining terms are 
at least linear in the spin parameter and of ©(C /i a/M), 
since they are all proportional to derivatives of which 
in turn is proportional to a/M [see Eq. ([M]) ]. 

The structure of Eq. (|5T|) thus reveals that if one were 
to employ the slow-rotation ap pro ximation, as in the 
computation of the MBH metric [22] (see Sec. lHI A]) , the 
leading-order evolution equation for the metric perturba- 
tion becomes 

C{De^) +2i?V>^i/)} = ^(V,V,^(i^o))*i?^^^^^ 

(54) 

Once this equation is solved for /i^ll^', one could use it to 
solve the evolution equation for which we have not 

presented here. In this way, one cou ld p roceed with the 
boot-strapping method described in [22] and in Sec. Ill Bl 
to build higher-order perturbations in the metric and the 
CS scalar field. 

In summary, GW emission in DCSMG is determined 
by the following equations. To first order, the metric 
perturbation is generated by the SCO through the GR 
relation in Eq. ([^5]) . leading to = 0(/i). To sec- 

ond order, the metric perturbation contains two contri- 
butions: a term h^l^f'' = 0{fj?) and a CS modification 
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to the first-order term /i^V'' = 0{ii(). The latter can be 
obtained by solving Eq. which requires knowledge 

of the leading-order perturbation to the CS scalar due to 
the SCO, which in turn is determined by Eq. 

once one has solved for through Eq. ([49]) . In this 



paper we shall only employ the first-order terms, h^^C^\ 
to model GWs, and thus modifications to the waveforms 
arise exclusively due to corrections to the trajectories. 

A proper treatment of the EMRI problem would also 
require the use of the first-order perturbations to es- 
timate backreaction, self-force effects on the SCO tra- 
jectory (see e.g. [S^l), which involves calculations that 
have not been carried out in full generality even in 
pure GR. Only through the adiabatic approximation has 
there been studies of such back-reaction for spinning 

MBHs m im, s m, m, s m, m m, p, where 

one assumes that only GW dissipation is important. 
The purposes of this paper, however, is to elucidate the 
main differences between GR and DCSMG dynamics for 
IMRI/EMRI systems, and thus, calculations with back- 
reaction are far too complex. Instead, we will resort to 
the so-called semi-relativistic approximation, which we 
describe in the next subsection. 



B. The Semi-Relativistic Approximation 

First introduced by Rufhni and Sasaki [50| . this ap- 
proximation assumes that the SCO moves along geodesies 
of the MBH background geometry and, more impor- 
tantly, that the SCO emits GWs as if it were moving 
in flat space. These are two strong assumptions but 
they give a good qualitative picture of the dynamics. 
In fact, it has been recently shown [s^ that in combi- 
nation with some post-Newtonian prescriptions for ra- 
diation reaction, semi-relativistic like calculations lead 
to waveforms with high overlaps with Teukolsky-based 
waveforms (waveforms were the GW fluxes for the adi- 
abatic approximation are calculated from solutions to 
the Teukolsky equation for perturbations of a Kerr back- 
ground). The geodesic motion assumption is justified for 
EMRIs because the more extreme the mass ratio, the 
closer one is to the point-particle limit, and hence the 
closer the motion is to geodesic. The flat-space assump- 
tion allows us to neglect the curvature of the background 
in the generation and propagation of GWs from the SCO 
to the observer in the radiation zone. 

Following the work of Ruffini and Sasaki [s^l , the im- 
plementation of the semi-relativistic approximation is as 
follows. Once the geodesic equations have been solved 
for the orbital trajectories, the metric perturbation is ob- 
tained from Eq. (|49p , but assuming that the background 
metric g^^ is the flat spacetime metric. In this way, the 
solution to the perturbative equation is given b y th e grav- 
itational Lienard-Wiechert potentials (see e.g. [8^[83| for 
the derivation of the electromagnetic version of these po- 
tentials), which contain all the information necessary to 
reconstruct the gravitational field generated by the SCO. 



In this paper, we shall employ a formulation of the 
semi-relativistic approximation in which the first assump- 
tion (geodesic motion with respect to the MBH back- 
ground geometry) is still employed, but the second one is 
slightly modified. More precisely, we shall still assume 
that GWs are emitted as if in flat space but, follow- 
|90|, we solve the perturbative equations using a 
multipolar expansion, which is determined by the source 
multipole moments and which is valid for objects with 
arbitrarily strong internal gravity (see [Slj for a detailed 
discussion of these expansions) . 

The multipolar treatment of gravitational radiation is 
truly a slow-motion approximation as the series is trun- 
cated at a finite multipole order. For a compact-body 
binary system, the error scales with a high-power of the 
orbital velocity, which for EMRIs can typically be larger 
than half the speed of light, and may lead to the loss of 
certain relativistic features in the waveforms (see eg. [9i| 
in the context of post-Newtonian theory and EMRIs). 
In spite of these drawbacks, the semi-relativistic approx- 
imation has been shown to capture the correct qualitative 
behavior of EMRI orbits and waveforms, and thus, it will 
suffice to present the main features of the CS modifica- 
tions to GW generation. 

Another possibility to estimate the GW emission is 
to use the weak- gravity/fast- motion formula derived by 
Press [92], whose implementation requires the same in- 
formation as the multipolar formula truncated after the 
mass octopole and current quadrupole. This formula has 
been used in (ssj . where it has been shown to provide 
accurate results as compared to Teukolsky-based wave- 
forms. We shall not employ this formalism here, however, 
choosing to work instead with the standard multipolar 
decomposition. 

We have shown that up to order 0{fi() the perturba- 
tive GW generation equations for IMRI/EMRI systems 
have the same form as in GR. Any difference in the wave- 
forms will then arise due to differences in the geodesic 
equations due to the modified Kerr background. In 
the transverse-traceless (TT) gauge, defined by Eqs. ([9]) 
and pop and in which the dynamical degrees of freedom 
have been isolated, the metric perturbations describing 
the GWs emitted by an isolated system are given by [51| 



^ £lr ^^JA,_J.'tret) ^A, 



i=2 ^ ' 

where we have adopted the multi- index notation of [Slj: 
is a multi-index (sequence of £ indices: ai ■ ■ ■ a^) so 
that Na — Ua ■ ■ ■ n,, is the product of i unit vectors 

= x^/r that point in the direction from the source 
toward the observer, with r = {SijX^x^Y/'^ is the flat- 
space distance from the source to the observer; there is 
summation over repeated indices independently of their 

ijk 



location; £,„ j, is the 3-dimensional flat-space Levi-Civita 
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tensor. 

The mass and current multipole moments of the 
source, and J4 , are given by [sij 



(56) 



where = x 



■ x^ and p is the mass density of the 
source, which in our case corresponds to the mass density 
of the SCO, given by 



p{t,x') ^m8'^^'>[x' - z\t)], 



(57) 



where S^'^) denotes the Dirac-deha distribution, z^{t) is 
the spatial trajectory of the SCO, and = dz^/dt its 
spatial velocity with respect to the coordinate time t. In 
Eq. (j55p the mass and current multipole moments have 



to be evaluated at the retarded time L 



t — r and the 



superscript {£) denotes their ^-th time derivative 

The projectors stf in Eqs. ([55)) and tt in Eq. ([55)1 are 
symmetric/trace- free and transverse/traceless operators. 
The latter is obtained by means of the projector orthog- 



onal to n% namely 



in the following way 



[A 



— P'^P'- A 



2 



A, 



(58) 



for the study of extreme-mass ratio gravitational-wave 
bursts. A different choice is to identify Boyer-Lindquist 
coordinates (r, 0, 0) with flat-space spherical coordinates 
and then introduce Cartesian coordinates through the 
familiar relations 



r sm f cos ( 



y = r sm U sm ( 



(61) 



This was the choice made in |55'| for the construction 
of EMRI kludge waveforms. Numerical experiments ^Q^l 
show that indeed different choices of Cartesian coordi- 
nates produce different waveforms. However, differences 
are significant only for relativistic orbits, and for those, 
these differences appear dominantly in the amplitudes 
but not in the phase. In this paper we shall employ the 
second choice of Cartesian coordinates for the computa- 
tion of the gravitational waveforms. 



V. GRAVITATIONAL WAVES FROM 
IMRI/EMRI SYSTEMS IN DCSMG 

In this section we study numerically the evolution 
of IMRI/EMRIs and their GW emission in the semi- 
relativistic approximation and in the context of the DC- 
SMG theory. We describe the numerical implementation 
and apply it to a number of test systems, presenting the 
trajectories and the waveforms associated with them. 



An orthonormal triad, with respect to the 

spatial flat-space metric and associated with n% can al- 
ways be constructed to build the polarization tensors via 



4j = PiPj - QzQi and e„- = 2pf-q- 



With these tensor projectors and Eq. we can ob- 
tain the two independent GW polarizations as follows: 



or equivalently 



~ 2 *i ■ 



(59) 



(60) 



In this paper we use this multipolar decomposition to 
compute the GW emission including terms up to the mass 
octopole and current quadrupolc moment. 

Another important question for the implementation of 
this version of the semi-relativistic approximation is the 
choice of coordinates. More specifically, the equations of 
motion (geodesies) are written in Boyer-Lindquist coor- 
dinates (see Sec. IIIIBp . whereas the formulae given in 
this section require Cartesian coordinates. Thus, one 
must decide how to construct Cartesian coordinates that 
cover the motion of the SCO and at the same time the 
observer in the radiation zone. One possibility would be 
to rewrite the background MBH metric and the equa- 
tions of motion in Cartesian-like coordinates such that 
at infinity the metric becomes the flat spacetime metric 
in Cartesian coordinates. This choice was made in 19011 



A. Numerical Evolution and Test Systems 

In the semi-relativistic approximation the numerical 
calculations can be divided into two stages. First, the 
computation of the SCO trajectory around the MBH and, 
second, the computation of the gravitational waveforms 
given the trajectory. The starting point for the com- 
putation of the trajectory, namely [t{T),r{T),6{T),(l){T)] 
are the geodesic equations [Eqs. ([50)1 and (P?]) ]. Clearly, 
when the CS parameter ^ is small the dynamics will be 
close to GR. Since there are already bounds on this pa- 
rameter from binary pulsar tests [22], albeit weak, we 
shall be here interested only in small CS deformation of 
GR EMRI/IMRI waveforms. 

As in GR, for bound orbits the equations for r(r) and 
9{t) present turning points which are difficult to treat 
numerically. This can be avoided by introducing the fol- 
lowing alternative variables (see e.g. jlHl for details): 



pM 



1 + e cos -0 



cos^ = cos2(6'„,J cos2(x) , (62) 



where p is the semi-latus rectum^ e is the eccentricity, and 
as we mention before 9 G (6'„i„,7i' — 6'„,i„)- The geodesic 
equations then become evolution equations for the vari- 
ables {ip,Xj4') with respect to Boyer-Lindquist time i, 
where one uses the dt/dr equation [Eq. ()30p ] to transform 
from proper to coordinate time. The form of the ordi- 
nary differential equations (ODEs) for [■0(i), x(t), (/)(i)] in 
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GR can be found from the developments presented, for 
instance in [s^ . 

From a numerical standpoint, ODE solvers for 
[ipit), x(^)j 4'i't)] require knowledge of the turning points, 
which will be here CS modified. In particular, the radial 
turning points arc modified, since the ODE for ip{t) is 
corrected by new CS terms that arise in the dr/dr equa- 
tion. In contrast, the polar turning points are not modi- 
fied, since these arise from the dO/dr equation, which is 
not CS corrected. All of this suggests that prescribing a 
set of initial constants of motion (i?, L, Q) leads to three 
orbital parameters (p, e, t), which are not identical to the 
ones that one would obtain in GR. 

Let us define these orbital parameters more carefully, 
since they play an important role in our calculations. The 
semi-latus rectum and the eccentricity parameter, p and 
e, are defined such that the apocenter ^ and pericenter 
''pori (two of the turning points of the equation for dr/dr) 
have the familiar Newtonian expression 



pM 
1-e 



pM 

1 + e 



(63) 



The third orbital parameter, the orbital inclination with 
respect to the equatorial plane i, is defined as tan l = 
Q/L, but it is sometimes more convenient to describe it 
in terms of the angle 9^^^ = 7r/2 — sign(L)6'j^;^, which is 
directly related to the turning points of 9{t). Clearly, the 
fimctional form of [p, e, 9^^^^) in terms of (_£, L, Q) differs 
in GR and in DCSMG, and thus, if one prescribes a set 
of constants of motion (Eq, Lq, Qq) in GR the resultant 
orbit will not possess the same orbital parameter as its 
DCSMG counterpart. In the results presented below, we 
compare trajectories that have the same orbital param- 
eters, (p, e,0;^J, instead of comparing trajectories with 
the same constants of motion {E, L, Q). 

Regarding the numerical integration of the ODEs for 
\tlj{t), x(t), Mt)], we use the Bulirsh-Stoer extrapolation 
method [94] as the evolution algorithm (details on this 
method and its implementation can be found in oH, |9^). 
Our numerical code is complemented with a number of 
routines to construct Boyer-Lindquist and Cartesian-like 
coordinates and their associated time derivatives. 

Waveforms are computed with the multipolar expan- 
sion of Sec. IIVBI up to the mass octupole and current 
quadrupole, which requires as input the trajectory a:'(i), 
the spatial velocity v'^{t), the spatial acceleration a'(i), 
and its time derivative, the jerk, f{t) = da'^{t)/dt. These 
time derivatives are computed using a finite differences 
differentiation rule with nine points. The error in the 
computation of the jerk scales with the time step. At, as 
(At)*, well within the accuracy range that we need for 
our calculations. 

We have here evolved the geodesic equations for a total 
time of T = 5 X lO'^M, obtaining on the order of 10^-10'' 
cycles. The type of geodesic is prescribed in terms of 
the following orbital parameters: the pericenter distance 
r ^^j, the eccentricity e, and the inclination angle 9^_^^. We 
implemented a double-bisection algorithm to guarantee 



that the turning points in both the GR evolutions and 
the DCSMG ones correspond to orbits with the same 
orbital parameters (j),e,9^^J, up to an accuracy of one 
part in 10^''. 

The MBH background metric of Eq. ^ is fully deter- 
mined by the choice of the mass M and spin parameter 
a = J/M, where J is the spin angular momentum, plus 
the ratio ^. The SCO is characterized by its mass m, 
which then also defines the mass ratio fi = m/M . As we 
have mentioned, the type of orbit is fully determined by 
the choice of orbital parameters {T^^^^,e,9^_^^), while we 
still need to fully determine the initial conditions of the 

motion by prescribing (i/'o, Xo> ^^o) = (V'(io)> xC^o): <^(*o))> 
where is the initial time for the numerical evolution of 
our system of ODEs. We shall here always set (3 = a, 
M = M, where A/, — 4.5 x lO^Af© is approximately the 
mass of the presumable MBH at the center of the Milky 
Way ■ For the mass of the SCO we choose m = 35 A/0, 
which leads to the following mass ratio: /i ~ 7.8 x 10~^. 

The particular test systems that we consider in this 
paper are defined as follows: 

(A) : a — 0.1 M and ^ = O.lAf^, with orbital parameters 

(V..e,a„J = (12Af,0.2,0.1), 

(B) : a ~ 0.2 M and ^ = 0.2Af'*, with orbital parameters 

(V.i'e,0.„J = (8A/,O.4,O.2), 

(C) : a = 0AM and ^ = 0.4Af'^, with orbital parameters 

(V.i'e,0.„J = (6Af,O.6,O.3), 

while the waveforms are measured by observers located 
on the z-axis at a distance of approximately 8 kpc 
(roughly the distance from the Solar System to the cen- 
ter of our galaxy). Clearly, the test systems have been 
ordered from the least relativistic to the most relativis- 
tic. We are forced to choose a and ( such that the 
slow- rotation approximation a/M 1 and the small- 
coupling approximation C ^ 1 used for constructing the 
MBH background '^2\ are satisfied. In this sense, the er- 
ror in the MBH background metric scales as C{a^/M^), 
which for systems (A), (B), and (C) implies errors in the 
MBH background of 0(10"^), 0(10-^) and 0(10-^) re- 
spectively, relative to the true DCSMG waveforms for an 
exact, spinning MBH background. 



B. Orbital Trajectories 

The test systems described above present rather simi- 
lar orbital behavior. Generically, there is a stage of zoom- 
whirl, where the particle spends several cycles close to the 
pericenter radius, r^^^^, followed by a stage where it or- 
bits at a larger radius r, close to r^^^. Moreover, there is 
generically both in-plane and out-of-the-plane precession. 
The orbit produced by test system {B) in GR is shown in 
Fig. [21 where we have plotted only the last time interval 
of duration 17500Af of the geodesic evolution. From this 
figure we can observe the generic Lense-Thirring preces- 
sion as well as the two stages described above. 
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FIG. 2: Three-dimensional depiction of the last cycles of sys- 
tem B before the evolution was stopped. 



The CS correction to the background has a clear ef- 
fect on the trajectories of the bodies. Figure [3] plots the 
projection of the orbital trajectories for system C onto 
the x-y plane for the CS background (black line) and the 
Kerr background (light gray line) and the last 17500M 
of the geodesic evolution. Notice that the orbital tra- 
jectories have dephased significantly by the end of the 
evolution. 



20 - 




FIG. 3: Two-dimensional projection onto the x-y plane of 
system (C)'s orbital trajectories in a Kerr background (light 
gray) and in the CS corrected Kerr solution (black). 



From this figure one cannot discern whether the CS 
trajectories trail or anticipate the GR ones, but this 
can be assessed by plotting the difference between the 
GR and CS evolved angles Xi that determine com- 
pletely the geodesic evolution. These angles are the 
Boyer-Lindquist azimuthal angle plus the two angles ^ 
and X associated with the other Boyer-Lindquist coordi- 



nates r and 9 through the relations in Eq. ((62)) . Figure |4] 
plots the dephasing in tp, where we observe that, in all 
panels, the difference Sip := ipcs — "^gr is initially close 
to zero and then it decreases linearly on average (over 
a number of cycles). Also we observe that, as expected, 
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FIG. 4: Dephasing of the if? coordinate as a function of time. 
The top panel corresponds to system (A), the middle one to 
system (B) and the bottom one to system (C). 



system (C) presents the biggest dephasing, followed by 
system (B) and then by system (A) (note the change of 
scale in the y-axis). Similar linear dephasing trends are 
observed in Sx and 6 (p. 

For weakly gravitating systems, we find that Sip is dom- 
inant and negative, while Sx —S(p, with the latter pos- 
itive. This implies that the CS orbits of systems (A) and 
(B) are overall trailing the GR solution. This scenario is 
a bit more complicated for system (C), where Sip <^0 is 
still dominant, but both Sep and Sx are positive, implying 
that the CS orbits still trail the GR one in radius, but 
anticipate it in angles and cp. 



C. Gravitational Waveform 

Orbital trajectories can give us a sense of the CS effect 
on geodesies, but the true observables are gravitational 
waveforms. Figure [5] plots the difference in the wave- 
forms computed with CS trajectories and GR trajecto- 
ries. For reference, the maximum magnitude of the GR 
GW polarization h_^^ is approximately \h^\ < 8 x 10^^^ 
[system {A)], < 1.25 x 10'^^ [system {B)], and 

\h_^_\ < 2.25 X 10'^^ [system (C)]. At the end of the 
simulation (~ 128 days of data using M = M,), we ob- 
serve that system (A) has dephased by as much as 0.3% , 
while system (B) has dephased by 16%, and system (C) 
by 90% relative to the maximum amplitude of the respec- 
tive GR waveforms. Similar behavior is observed for the 
other polarization. 

Another measure one can construct to observe this de- 
phasing is the time-dependent overlap Q{gr, cs), defined 
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FIG. 5; Dephasing of the GW polarization as a function 
of time. The top panel corresponds to system (A), the middle 
one to system (B) and the bottom one to system (C). 



as 



0{A,B):=hih^ + hih^, (64) 

which we can normalize by using the quantity 

Ojv (gr, cs) := a/ 0(gr, gr) O (cs, cs) . (65) 

This measure is related to the integrand of the overlap 
commonly used in gravitational wave data analysis. Fig- 
ure[S]plots Q{gvi, cs)/Qpf as a function of time. Observe 



0.999995 




wave astronomy can be when considering geodesic mo- 
tion about a spinning MBH, since even a rather strong 
curvature correction to the action leads to gravitational 
waves that are virtually indistinguishable from GR ones, 
unless the SCO is close to the light-ring or the coupling 
constant of the curvature correction is quite large. 

By considering system (C), the most relativistic one, 
one can ask how small can C be in order to cause a loss 
in the normalized dephasing measure of about 5% after 
10^ M of evolution [Eqs. ([M]) and ([HS])]. The answer can 
be found in Figure [71 which shows the average of the 
normalized dephasing measure as a function of time for 
system (C), but for the following values of the CS pa- 
rameter C: C = (solid black), ^ = 0.2M^ (dotted 
blue), C = O.IM^ (dashed red), ^ = O.ObM'^ (dot-dashed 
green), and ^ = 0.025Af'* (dot-dot-dashed violet). Ob- 
serve that over 4 months of data, the CS correction could 
lead to a significant dephasing even for a CS parameter 
of 0(10^^). On the other hand, if one integrates inco- 
herently over 3 week segments of the data, then the CS 
correction would lead to a dephasing only if the CS pa- 
rameter if of 0(10"^). In this case, however, one would 
have trouble connecting three-week segments together if 
the CS correction is not taken into account. 
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FIG. 7: Normalized dephasing measure 0{gr, cs) as a func- 
tion of time for system (C) using: ^ = (solid black), 
^ = 0.2M'' (dotted blue), ^ = 0AM* (dashed red), ^ = 
0.05M* (dot-dashed green), and ^ = 0.025M'' (dot-dot- 
dashed violet). 
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FIG. 6: Running average of the normalized dephasing mea- 
sure 0(c:f^)Cs) over twenty consecutive time steps as a func- 
tion of time. The top panel corresponds to system A, the 
middle one to system B and the bottom one to system C. 



that after only three weeks of data (roughly 10^ Af of evo- 
lution) , system {A) and {B) will probably not be able to 
distinguish between GR and CS geodesies. On the other 
hand, system (C) will probably be able to distinguish this 
difference, because it experiences a deeper gravitational 
potential. This plot shows how robust LISA gravitational 



Another important question is how this bound com- 
pares with bounds obtained from binary pulsar tests. 
The precession of the perigee in binary pulsars can be 
used to argue that (^^^^ < 10^ km, as was done in Paper 
I. For inspirals, a GW test of DCSMG can be expressed 
as 



^^/^ < 6 X 10^ S^/"^ {M/M,) km , 



(66) 



where S here represents the accuracy to which we can 
measure C. This accuracy depends not only on the in- 
tegration time, the signal-to-noise ratio and the type of 
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orbit considered, but it should also be enhanced if one ac- 
counts for radiation-reaction effects, since they are also 
CS corrected (see Sec. |Vl|. Based on the results of this 
section, it is not ludicrous to expect that S ~ 10~^ or 
better depending on the system, while one also sees that 
IMRIs are favored over EMRIs due to the M'^ depen- 
dence. Thus, if one considers an IMRI with total mass 
M = 10"^ M0, and if one assumes S — 10^^, then the 
constraint can be better than the binary pulsar one by 
at least two orders of magnitude. Moreover, as we shall 
see in Sec. I VII radiation-reaction effects depend on the 
oscillatory sector of the CS scalar, which in turn affects 
GWs directly. Such an oscillatory sector of the theory is 
simply untestable by binary pulsar experiments. 

The results and arguments presented above suggest 
that GW observations with LISA could allow for inter- 
esting tests of DCSMG, but of course, these arguments 
should be backed up by a deeper and more thorough data 
analysis study. In this sense, the figures presented in this 
section should only be taken as an indication of the pos- 
sible dephasing and loss of overlap that one could expe- 
rience if CS theory is in play. A detailed study of the 
dependence of such a test on the spectral noise density 
curve, the location of the detector in the sky via the beam 
pattern functions, and the distance to the source, is thus 
of key importance, but postponed to future work. 

VI. GRAVITATIONAL WAVE PROPAGATION 
AND RADIATION REACTION 

The semi-relativistic approximation provides a sensi- 
ble approximation to the emission of GWs generated 
by the the motion of SCOs around MBHs, but so far 
radiation-reaction has been completed neglected. These 
effects can be incor por ated via the adiabatic approxima- 
tion [il. lisl. lirl. Ilsl.^?! . that is, the assumption that the 
changes in the constants of motion {E, L, Q) evolve in 
a time scale much larger than the orbital ones. In this 
scheme, one then argues that the system can be evolved 
using a geodesic evolution for a certain number of orbital 
cycles, after which the constants of motion are corrected 
by estimates that use balance laws (eg. the GR energy 
balance law provides a prescription of how to modify the 
orbital energy due to GW energy emission out to infinity 
and into the horizon). 

In order to obtain the balance laws one must under- 



stand the propagation of GWs in the underlaying theory 
and their associated, effective stress-energy tensor. This 
can be done in the framework of the shortwave approxi- 
mation, where one decomposes the geometry into a back- 
ground and an oscillatory part corresponding to GWs. 
Such a decomposition holds when the GW wavelength 
is much smaller than the typical length scale associated 
with the background curvature. 

In this section, we develop the shortwave approxima- 
tion for DCSMG and compute the effective stress-energy 
tensor of GWs in this theory. Based on this, and us- 
ing the symmetries of the background geometry (time- 
like and axial Killing vectors), we can establish balance 
laws that can allow us to introduce radiation-reaction ef- 
fects, in the adiabatic approximation, in the energy and 
angular momentum. For the case of the Carter constant 
the situation is more difficult as it is not associated with 
a Killing vector symmetry, but perhaps a two-timescale 
expansion could be performed, similar to that introduced 




A. The Short- Wavelength Approximation 

The short- wavelength approximation (SWA) was first 
introduced by Isaacson [9^ ilOQ] (see also [63 and ref- 
erences therein) to study the propagation of GWs be- 
yond the linear approximation. In the SWA one assumes 
that the typical GW wavelength, A^^, is much smaller 
than the curvature length-scale, TZ, such that A^^ <C TZ, 
where TZ~^ is the typical scale of the components of the 
Riemann tensor. In this way, we can split the spacetime 
metric into a background metric, g^^,, with associated 
curvature scale TZ, plus an oscillating metric perturba- 
tion describing the GWs, /i^^, with an associated length 
scale Aq„ [see eg. Eq. ©] and an associated amplitude 
/lew In this way, d^g^^ 1/7?. and dah^^ ^ ^Iqw/'^gw 

In practice, applying the SWA to any tensor implies 
separating its background part from its oscillatory part, 
which can be done by averaging over several wavelengths 
(•^Gw)- Such separation can be ac comp lished via the 
Brill-Hartle averaging procedure (99l . llOOj ]. which applies 
to convex regions of the background spacetime (regions 
where any two points can be joined by a unique geodesic) . 
Then, the average of the components of a given tensor 
T^i ' bi... is defined as: 



, rpai ■ 



>(:.):= / dVV^i(^/(x^a;)g^}(x,a;0•••g£(x,xO•••T;'C(x') 



w 



(67) 



where the integration takes place over a small region, W, sor of parallel displacement in the background spacetime; 
containing several GW wavelengths; g""^ a[ is the biten- and f{x', x) is a weighting function that smoothly decays 
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to zero as x and x' are separated by many wavelengths. 
The main properties of this averaging rule (for details 
see [i^) are: 

1. Background covariant derivatives commute [frac- 
tional errors are of ©(A^^/T?.)^]. 

2. Total divergences average out to zero [fractional er- 
rors are of 0{X^„/TZ)]. 

3. As a corollary of (1) and (2), one can freely inte- 
grate by parts. 

These rules shall be used heavily in the next section to 
obtain an expressions for the GW stress-energy tensor in 
CS modified gravity. 



-^l*R^'^R[h] + jh*RR 



(70) 



The detailed structure of Eq. ([55)1 as well the form of 
6TP[^,d] can be found in Appendix [X2] [Eq. (|MT|) ]. 
The object ''^'S'^a in Eq. (|70p denotes the perturbation 
of the Christoffel symbols and its expression is given in 
Eq. (1X91) . 

Using the averaging procedure described above, the 
field equations tell us how the nonlinear contribution 
of the waves shapes the background. This contribution 
can be written as an effective stress-energy tensor of the 
GWs, T^™ (the Isaacson tensor), allowing us to write 



B. Effective Stress-Energy Tensor for GWs 

In DCSMG, not only does the metric tensor oscil- 
late, but also the CS scalar field, sourced exclusively by 
the spacetime curvature [see Eq. ([5])]. When GWs are 
present, these will then generically induce oscillations in 
the CS scalar field, which forces us to also decompose 
into a background part, "(9, and an oscillatory part, z?, the 
latter induced by GWs. 

Let us now study the field equations in the SWA, which 
requires expansions up to second order in the oscillatory 
parts (see Appendix lA II for the basic formulae used for 
these expansions). The structure of the resulting equa- 
tions consists of oscillatory terms, linear in the ampli- 
tude, and coarse-grained terms that describe how the 
background is modified by GWs. This modification is 
produced by the averaged second-order term in the oscil- 
latory part, which refiects the nonlinear character of this 
effect. 

Then we set to zero the linear part in the wave am- 
plitude, /i(3wi obtaining in this way the equations that 
describe the propagation of the waves in the background. 
In GR, these equations are simply 
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(68) 
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(1) T^mat 



where we have used the Lorenz gauge [Eq.Q] and where, 
here and in the rest of this section, the superscript pre- 
ceding a given quantity denotes the perturbative order 
of this quantity with respect to metric perturbations (see 
Appendix [X] for more details). 

In DCSMG the propagation equations have the follow- 
ing form 



_<i)r(^)[^ 



" 2 S/ii/ ( 



^(<^'C^,[z9,/i] + C^,[^, 
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2k 



(69) 



(71) 



In GR, the Isaacson tensor is given by (see e.g. 



(72) 

which, using the expressions in Appendix |^ and work- 
ing in the Lorenz gauge [Eq. ([5])] supplemented with the 
traceless condition [Eq. (fTO|l ]. reduces to: 



(73) 



One can see, using the first-order equations and the aver- 
age rule, that this tensor is traceless and divergence- free 
within the errors produced by the approximations made 
using the SWA scheme: 



Tf^VrpGW n rpQVj\l' n 

6 ^ fill -'- 



(74) 



In DCSMG, the effective stress-energy tensor of GWs 
is given by 

= -2k |< <^'i?^,[/»] > -ig^, < ^^^R[h] > 

+ ^ (< '''C^A^,h] > + < ''^C^A^,h] >)}(75) 

Again, using the formulae in Appendix 13 restricted to 
the TT gauge of Eqs. ([9|) and (fTO|) . and the properties of 
the averaging (in particular the integration by parts) , we 
find that the terms corresponding to the CS correction 
vanish exactly^ yielding the same expression as in GR (for 
a guided explanation of this fact, see Appendix [A|. That 
is, the effective stress-energy tensor of GWs in DCSMG 
(which enters in Eq. (|7ip ) is simply given by Eq. (j72[) 
or ([73|) , which implies that the backreaction in the orbital 
motion due to GW emission is essentially as in GR. 

We can now look at the equation for -d with the second- 
order corrections. This is a linear equation in and hence 
all the modifications come from second-order terms in the 
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GW. One the finds that 

pad = -^(^l-l{h-0h^f,)^ *RR 

- |e«/3M.('^'i?"''''"'^'i?p.^'') . (76) 

As we can see there are two second-order corrections to 
the background value of the CS scalar field: one that 
renormalizes the value of a, < h^^h^j^ > ~ ^^'^ 
another that scales as ah^^/X^^. 

C. Balance Laws and Radiation- Reaction Effects 

By means of the SWA we have studied the propagation 
of GWs and of oscillations induced in the CS scalar field. 
We have also seen how these oscillations affect the back- 
ground values of the metric tensor [Eq. ([TTjl ] and of the 
CS scalar field [Eq. ((76|) ]. Of particular importance is the 
equation for the modification of the background geome- 
try, which contains the effective stress-energy tensor of 
the GWs, which we have found has the same form as in 
GR. From the previous development it is clear that from 
Eq. (j7ip we have the following conservation equation 

V^T =0 T T""' -f r'^''H?91 -I- T'^" (77) 

where we have used that the background Einstein ten- 
sor is divergence- free and the equations of motion for the 
background scalar field d are satisfied. The new ingredi- 
ent with respect to GR is the appearance of the stress- 
energy tensor associated with the oscillations induced in 
the CS scalar field. 

Here, the background clearly corresponds to the MBH 
geometry and CS background scalar field is given in 
Eqs. ([23|) and ([24|) . Using the Killing symmetries of the 
background, described by the vector fields and ip'^ [see 
Sec. lIIIB] . the vector fields 

- -r^,t" , = T^.i^" , (78) 

describe the total fiuxes of energy and angular momen- 
tum and are, by virtue of Eq. (|77p and the Killing equa- 
tions, divergence-free vector fields. Therefore, we can 
obtain balance laws for the energy and angular momen- 
tum by considering a spacetime region V with boundary 
dV and integrating over it the divergence- free conditions: 
f'^l^ = and Jt'i^ = 0. 

For the study of IMRI/EMRIs one can take V = 
{{t,r,6,(j)) , such that < t < tj and r^j < r < r^}, 
where the inner radius, r^j^ can be taken to be close to 
the horizon location and the outer radius, r^, can be 
taken to be effectively close to infinity (or at least far 
away from the SCO). In this way dV is composed of two 
cylinders (one aX r — and the other one at r = r^) 
of finite size, limited by two slices, one at an initial time 
t — and the other one at a final time t ~ t^. 

The balance laws for energy and angular momentum 
tell us that the change in the constants of motion {E, L) 



is given by the amount energy/ angular momentum, of 
the GWs and of the CS scalar field, fiowing away from V 
to infinity [r = rj) and into the MBH horizon {r — r^). 
In this way, one can account for the radiation-reaction 
effects in the adiabatic approximation. Note that this 
only takes care of the changes in the energy and in the 
angular momentum of the SCO, while the question of 
how to correct the third constant of motion, the Carter 
constant Q, in a general way remains still open (see [H, 
m, m, El] for recent advances on that question) . 

The flux of energy (or energy luminosity) in GWs going 
towards infinity is given by 

- - fim / dn T^T , (79) 

dt r^oo I q2 

where = x^/r is the unit normal to the surfaces 
r — const, near spatial infinity, where the background 
geometry is essentially flat. Although the expression for 
the energy luminosity [Eq. (|79p ] is formally the same as 
in GR, the luminosities are actually different because the 
geodesic equations are indeed CS corrected by the mod- 
ifled Kerr background given in Eq. 

To illustrate this difference, let us consider the expres- 
sion 

^ = .[hl + hl), (80) 

where and hx are the GW plus- and cross- 

polarizations, which one can obtain from Eqs. (|73p 
and ([79]) . Figure [8] plots the difference between the en- 
ergy flux computed from the GW emission of a SCO mov- 
ing in geodesies of Kerr and a SCO moving in geodesies 
of the CS-modified Kerr metric [Eq. ([SS])]. for systems 
(A) (top panel), {B) (middle panel) and (C) (bottom 
panel) as a function of time. For reference, the maximum 
magnitude of the energy flux computed from geodesic 
waveforms about Kerr are 3 x 

10-37m"2 foj. system (A), 
3 X IQ-^^M-'^ for system (B) and 2 x IQ-^^AJ-'^ for sys- 
tem (C). Observe that, after approximately four months 
of evolution, the fractional correction to the luminosities 
by the CS correction is about 1% for system (A), 10% for 
system (B) and about 50% for system (C). As expected, 
the more relativistic the system, the larger the CS cor- 
rection to the waveforms, and therefore, the larger the 
effect on radiation-reaction. Such corrections could have 
a significant impact on kludge- type waveforms, when one 
takes into account that the GW energy lost must also 
be supplemented by the energy emitted by the CS scalar 
field [see the conservation equation, Eq. (fTZl) ]. To include 
this contribution, one would need to solve the evolution 
equation for the CS scalar field [Eq. ([5D)) ]. using as input 
the gravitational waveform computed without radiation- 
reaction. 

The machinery developed in this paper allows, for the 
first time, the inclusion of these radiation-reaction ef- 
fects in the orbital evolution of EMRIs and IMRIs in a 
realistic, alternative theory of gravity. With these de- 
velopments, one can now properly construct alternative- 
theory, kludge waveforms. 
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t/M 

FIG. 8: Energy flux difi'erence (in units of M'"^) between the 
flux computed with gravitational waves produced by geodesies 
about Kerr and about the CS modified Kerr background for 
systems {A) (top panel), (B) (middle panel) and (C) (bottom 
panel) as a function of time. 



VII. CONCLUSIONS 

We have studied the effects of dynamical CS modi- 
fied gravity on the orbit and GWs generated by IMRIs 
and EMRIs. The semi-relativistic approximation was 
employed, through which we modeled the SCO trajec- 
tory via the geodesic equations, while we treated GWs 
through a multipolar generation formalism. 

We began by showing that test particles follow 
geodesies in DCSMG. We then proceeded to explic- 
itly calculate the modified geodesic equations for a test 
particle in a generic orbit around a CS modified Kerr 
MBH 01 . We showed that the fundamental frequencies 
of the modified background are CS modified, as well as 
the relations between the MBH multipole moments and 
its mass and spin, although the latter incurs corrections 
at £ = 4 multipolar order, which might be small to be 
detected by future GW observations with LISA. 

We then continued with the study of GW generation 
in the modified theory. We saw that the GW generation 
formula based on the multipolar expansion is corrected 
to leading-order by second-order terms of ©(/j-C), which 
we neglect here. These GW generation formalism is also 
modified by GR radiation-reaction effects, which scale as 
the square of the mass ratio and which we also neglect 
here. 

We then evolved geodesies in this modified background 
and compared their associated GWs to those generated 
by geodesies on a pure Kerr background. Generically, we 
find that these waves dephase after a certain number of 
cycles that depends on the strength of the CS coupling 
parameter, as well as on how relativistic the geodesic 
under consideration is and the signal to noise ratio of 
the event. For small CS parameters, such a dephasing 
is small, leading to mismatches that will not affect GW 
detection. These mismatches, however, might affect GW 



characterization, leading to erroneous conclusions about 
the orbital parameters of the system. 

EMRIs and IMRIs can thus be used to probe DCSMG, 
both by constraining the CS coupling strength to un- 
precedented levels, as well as to sample the oscillatory 
dynamic sector of the CS scalar field. Moreover, the 
study presented can be viewed as an explicit example 
of the effect a non-Kerr background could have on GWs, 
where in this case the non-Kerrness arises from a well- 
defined and physically well-motivated alternative theory. 
Since this modified theory introduces a higher-order cur- 
vature correction, modifications to the Kerr metric are 
only dominant in the strong-field and cannot be sampled 
by the waveform's first multipoles. 

The analysis presented here employs some approxima- 
tions that leads to certain inaccuracies. These inaccura- 
cies are associated with: (i) neglecting radiation-reaction 
in the dynamics, which scale with the square of the mass 
ratio; (ii) neglecting CS modifications to the GW emis- 
sion, which scale with the product of the mass ratio and 
the CS coupling constants; (iii) neglecting CS modifica- 
tions to the MBH background metric that are of third- 
order in the spin parameter and the small-coupling ap- 
proximation. Whether one effect is dominant over the 
others depends on the strength of the CS coupling con- 
stants, which is at present unknown. 

Future work could concentrate on the inclusion of 
radiation-reaction in an adiabatic fashion to the wave- 
forms presented here. Such a task in the modified the- 
ory is rather straightforward since these dissipative ef- 
fects are CS corrected only through the emission of CS 
scalar field. One could then supplement the standard 
GR expressions for the rate of change of the geodesic 
constants of motion by these CS corrections to allow for 
a true, kludge-like, inspiral. Although the CS correction 
does not modify radiation-reaction through changes in 
the Isaacson tensor, dissipative effects will exentuate the 
observed dephasing. 

Once such radiation-reaction effects are included, a 
thorough data analysis study should be carried out to 
confirm or not the expectations presented in this pa- 
per. Of particular importance will be the inclusion of 
the beam pattern functions and the motion of the de- 
tector in the sky, which h ave been shown to be relevant 
for parameter estimation [lOll lioj . Il03l . [lO^ . Moreover, 
the computation of the proper overlap and Fisher matrix 
would include the spectral noise density curve of the de- 
tector, as well as details of the sensitivity of the modified 
waveforms on the parameters of the alternative theory. 

Another avenue for future research is the inclusion of 
first, leading-order CS correction to the GW emission 
framework. We have provided here a map, albeit in- 
volved, to construct such a correction which is based 
on black hole perturbation theory. Another possibil- 
ity would be to study GW generation in the context 
of a post-Newtonian /post-Minkow s kian exp ansion and 
asymptotic matching [loi, [lOi, [lol|Iol Eoll • Although 
this formalism is formally valid in the slow-motion ap- 
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proximation, so is the truncated multipolar decomposi- 
tion used here. 

Last but not least, an exact metric that describes spin- 
ning MBHs with arbitrary spin angular momentum in 
DCSMG is still lacking. The results presented here are in 
principle valid only in the slow-rotation approximation. 
Generalizing such results to higher order in a perturba- 
tive fashion could be dramatically difficult, since to next 
order the scalar stress-energy tensor will also contribute 
to the modified field equations. Perhaps, a full numerical 
study of the collapse of a scalar field with angular mo- 
mentum would be appropriate to disentangle the effects 
of DCSMG on rapidly rotating MBH backgrounds. 
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APPENDIX A: FORMULAE FOR THE SHORT 
WAVELENGTH APPROXIMATION 

In this Appendix we provide some of the formulae that 
we have computed and used to derive the form of the ef- 
fective stress-energy tensor of GWs for DCSMG, using 
the SWA. However, these are generic perturbative ex- 
pansions (induced only by perturbations of the metric 
tensors) of geometric quantities which can be used in 
other perturbative schemes, although the interpretation 
of the different terms would be different. 



1. Second-order expansion 

Given the splitting of the spacetime metric given in 
Eq. ([5]) , we need to perform an expansion in h^^^ to second 
order of the different geometric quantities involved in the 
CS modified gravitational field equations. The inverse of 



the metric is: 

= gA"' - h^"" + h^phP"" + 0{h^) , (Al) 

where 

The expansion of the metric determinant, g = det(g^j^), 
is then given by: 

g = gf^l + h+^{h^-hP-'h^,)+0{h^)^ . (A3) 

From here we can expand the completely antisymmetric 
Levi-Civita volume four-form: 

e = ^/^<5"l^3 (A4) 

to get 

(A5) 

The Christoffel symbols can be written exactly as follows: 

r;^, = f^p, + , (A6) 

where is 

This tensor can be expanded in /i^^ yielding 

5;, = '^>5;, + <='5;, + o(;i^), (A8) 

where 

The introduction of the tensor S"^^ is convenient since it 
can be used to simplify the expressions for the expansions 
of the Riemann tensor and derived geometric quantities. 
An exact expression relating the Riemann tensor of the 
spacetime to the background one through the tensor 
is: 

R^pfiu = pnu ~ 2 'S'^i^iiy] + 2 '5'"[p5'^]^ . (All) 
Therefore, an exact expression for the Ricci tensor is: 

— + S'piy\a " ^ap\iy + ^ap^pi, ~ Spj,Sp^ . (Al2) 

From these expressions of the Riemann and Ricci tensors 
we can immediately write down their first- and second- 
order terms in the perturbative expansion in /i^^. For 
the Riemann tensor they are: 

^^''R" ppv = ~2 '^'S'^j^Ij,] , (A13) 
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(='i?«,_,„„ = -2 + 2 <^>S';^[^ '^'S';;]^ , (AM) The Einstein and Cotton tensors have more compHcated 



and for the Ricci tensor: 



expressions. The first- and second-order expansion of the 
Einstein tensor are: 



(1) D _ (1) ca _ (1) ca ( 



(2) p _ (2) (2) i(i)ca (1) qP (1) ca (l) cP 



(A16) 



(A17) 



"pa 



"pa 



pa J r\ 



"pa 



pa) 



(A18) 



The expansion of the first piece of the Cotton tensor, C^^, that these are expansions only in the metric tensor, the 
is determined by the foUowing expressions [we remark CS scalar field is assumed here to be fixed]: 



(1)^1 _ a - aafi{{l)-n (i) qP D \ i ,9 (i), aafS n 



(A19) 



2)/-rl — ,9 J aal3((2)T^ (i) qP (i) r> (2) c-P r> \ i ,9 (D, aal3 ( (l) n (1) qP D \ 

^t^iy - ^\a^{p \ ^i^)l3\a~ '^u)a ^Pp ~ '^u)a^Pp j + ^\a ^(^1 \ '^i/)a^/3p/ 



I ,9 (2), aaf3 n 



(A20) 



where (^'e'^^"/' and (^'e'^^"^ are the first- and second- The expansion of the second piece of the Cotton tensor, 
order expansion terms of e"^ f/^^ . In particular, the first- C^^^, is given by: 
order term is given by 

<^'e/"'3 ^-/ij6-/"'3-H 2 e-/A'"/i''" + ^/i (A21) 



(1)^2 _ i ,9 

^pp — 2 \aT 



- CTa/3 (1) pr _L aa0T)T 



v)afi 



(1) QP e '^"Z' nr 



(A22) 



(2)(-<2 _ 



pu 



(TQ/3 (2) nr 



R^ v)aP + ' '^(p 



(i)c, aaP{l)j^T ^ a_ (2)^ cra/3 pr 



i.)a/3 ':(p 



= (P 



cra^ (1) j^r 



(1)^ aap-p^T 



u)o,f3 ^ '=(p 



1 



,9 (1) CP J aaf} T)T 



(A23) 



r 



2. Propagation Equations in the SWA 



The propagation of oscillations in DCSMG using the 
SWA approximation is described by the first-order equa- 
tions given in Eq. (l69|) . Here, we analyze the differ- 
ent terms that appear in this equation adopting the TT 



gauge defined by Eqs. ^ and (|T0|) . We also make use 
of the SWA to simplify the form of some of the terms. 
Then, taking into account that 



<^'^S^ = 0, g'^''<^)^;, = 0, (A24) 
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the first term in (l69l) can be written as 



-h : P + ( 



(1) D _ (1) CO 



From Eq. ([7]) we have 



(A25) 



(A26) 



(A27) 



where ^^-^^ denotes the amplitude of the oscillations in 
the CS scalar field. From Eqs. (|M9l) and (|J22l) we can 
write: 



(1)^2 



9 - aaf3 (1) CP 



(A28) 



(A29) 



where denotes the magnitude of the background value 
of the CS scalar field. The last term in Eq. ((69|) that we 
have to look at is the variation in the stress-energy tensor 
of the CS scalar field, 5T'j£}[d, which is given by 



^tW[z9,^]=2/3 



(A30) 



Then, the first-order propagation equations [Eqs. 
assuming for simplicity that the waves travel in the ab- 
sence of matter fields, can be written as 



'5° 



9 - crafi (1) CP 
^|(T^(p '^v)l3 



Q I ^GW '^b 



^Gw'^ ^^Gw'^ ^GW'^ 



(A31) 



As it is clear, this propagation equation for /i^^ is coupled 

to the propagation equation for t9 [Eq. fTO]) ]. Notice that 
the right-hand side of the equations can be written as a 
total derivative, both in GR and in DCSMG. Also note 
that the second term in the right-hand side is to leading- 
order of 0(/i^„^,/(A3^72)). 
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